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Abstract 

This work introduces, examines and compares several quaternion 

normalization algorithms, which are shown to be an effective stage in the 
application of the additive extended Kalman filter (EKF) to spacecraft attitude 
determination, which is based on vector measurements. Two new normalization 
schemes are introduced. They are compared with one another and with the known 
brute force normalization scheme, and their efficiency is examined. Simulated 
satellite data are used to demonstrate the performance of all three schemes. A 
fourth scheme is suggested for future research. 

Although the schemes were tested for spacecraft attitude determination, the 
conclusions are general and hold for attitude determination of any three 
dimensional body when based on vector measurements, and use an additive EKF for 
estimation, and the quaternion for specifying the attitude. 
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I. INTRODUCTION 


Attitude determination of spacecraft usually utilizes vector measurements 
such as Sun, center of Earth, star, and magnetic field direction to update the 
quaternion which determines the spacecraft orientation with respect to some 
reference coordinates in the three dimensional space [1,2,3]. These measurements 
are usually processed by an extended Kalman filter (EKF) which yields an 
estimate of the attitude quaternion [4-8]. 

Two EKF versions for quaternion estimation were presented in the 
literature; namely, the multiplicative EKF [4-6] and the additive EKF [5,7,8]. 
In the multiplicative EKF it is assumed that the error between the correct 
quaternion and Its a-priori estimate is, by itself, a quaternion that represents 
the rotation necessary to bring the attitude which corresponds to the a-priori 
estimate of the quaternion into coincidence with the correct attitude. The EKF 
basically estimates this quotient quaternion and then the updated quaternion 
estimate is obtained by the product of the a-priori quaternion estimate and the 
estimate of the difference quaternion. In the additive EKF it is assumed that 
the error between the a-priori quaternion estimate and the correct one is an 
algebraic difference between two four-tuple elements and thus the EKF is set to 
estimate this difference. The updated quaternion is then computed by adding the 
estimate of the difference to the a-priori quaternion estimate. 

If the quaternion estimate converges to the correct quaternion, then, 
naturally, the quaternion estimate has unity norm. This fact was utilized in the 
past to obtain superior filter performance by applying normalization to the 
filter measurement update of the quaternion [7]. It was observed for the 
additive EKF that when the attitude changed very slowly between measurements, 
normalization merely resulted in a faster convergence [7,8]; however, when the 
attitude changed considerably between measurements, without filter tuning or 
normalization, the quaternion estimate diverged. However, when the quaternion 
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estimate was normalized, the estimate converged faster and to a lower error than 
with tuning only. 

In the next section we introduce the additive EKF for attitude 
determination. The role of quaternion normalization in the additive EKF is 
explained in Section III. In Section IV we discuss the brute force (BF) 
normalization scheme and examine its performance. In the following sections we 
introduce the quaternion pseudo-measurement (QPM), and the magnitude 
pseudo-measurement (MPM). Test results of the application of all normalization 
algorithms discussed in this work to simulated Earth Radiation Budget Satellite 
(ERBS) data is presented in Section VII. In Section VIII we introduce the 
linearized orthogonal ized matrix (LOM) normalization scheme as a suggestion for 
future investigation. Finally, the conclusions of this work are discussed in 
Section IX. 


II. THE ADDITIVE EKF FOR QUATERNION ESTIMATION 

Attitude determination from vector observations using the additive EKF is 

explained as follows. Suppose that a sequence v^ m ^ 1=0, 1,2,... of vector 

measurements performed in body, b, coordinates are given. Given are also these 

vectors in the reference coordinate system r. Denote the latter vectors by v^ ^ 

i =0,1,2, The vector v L . is a column matrix whose elements are the 

— bm, 1 

components of a vector v measured at time t^ and coordinatized in the body 

coordinate system. Similarly, the corresponding v , vector is a column matrix 

r , i 

whose elements are the components of the same vector v coordinatized in the 
reference coordinate system. Our aim is to estimate the quaternion q which 
expresses the body attitude with respect to the reference coordinate system. To 
meet this end we define an effective measurement y as follows 


zi 


v, , - A(q)v 
-bm, i — -r, i 


( 2 . 1 ) 
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where A is the direction cosine matrix (DCM) which transforms vectors from r to 

A 

b, and where q is the latest estimate of q. The vector v, . which is a result 
- — -bm, i 

of a measurement, contains all the error associated with the instrumentation, 
such as instrument misalignments, scale factor error, bias, white noise etc. The 
vector v is taken from the almanac and is assumed to be perfectly known. We 
observe that when the measurement is error free and when the quaternion estimate 
is accurate, is zero. On the other hand, when these assumptions do not hold, 
then y^ is a, generally non-linear, function of the instrument and attitude 
errors. 

The measured vector v^ m ^ can be expressed as follows 


dv 

v = v + 

— bm, i -b, i de 


( Se + ) 


r , 

-b, 1 


( 2 . 2 ) 


where v^ ^ is the error-free value of v when coordinatized in the b system, the 
Jacobian matrix 


a 5v 
H . A “ 


e, l 


de 


b, i 


(2.3) 


is the sensitivity matrix of the error associated with the measurement v as 

-bm, i 

a function of the instrument errors. The latter are expressed as a sum of a 
narrow spectrum error vector, 5e, and a wide spectrum error vector n^, which is 
modeled as a white noise error vector. The vector e contains all the 
instrumentation errors mentioned before, while 5e denotes the difference between 
e and its compensation value which is the latest estimate of e denoted by e. 
Define 5q as follows 


= q - q 


(2.4) 


then 


A(q) = A(q+5q) 


(2.5) 


406 



therefore, based on the assumption that 5q is small such that q is close 
to q, A(q) can be approximated as follows 


Define 


Since 


4 9A(q) 

A(q) = A(q) + [ 

J=1 J 


5q, 


consequently 


4 3A(q) 

A ( q)v r>1 - A(q)v rtl - [ g— - 


J=1 


«q 


3A(q) 


G J ‘ G J - ' 


A Cq) - 


2 2 2 2 
q r q 2 _q 3 +q 4 


2( q i q 2 c ?3 < ?4^ 


2( q i q 3 + q 2 q 4 ) 


2(q l q 2 +q 3 q 4 ) 

2 2 2.2 

- q l +q 2" q 3 +q 4 

2(q 2 q 3‘ q l q 4 ) 


2(q l q 3 _q 2 q 4 ) 


2( q 2 q 3 +q i q 4 ) 

2 2 2 4 

- ,] 1-W q 4 


then 





A A A 


q l q 2 % 


- q 2 q l _q 4 

G ! =2 

q 2 - q ! q 4 

AAA 

(2. 10a) G 2 = 2 

q l q 2 q 3 

AAA 


q 3 " q 4 _q l 


q 4 q 3 " q 2 



AAA 


AAA 


“ q 3 q 4 q i 


q 4 q 3 -q 2 

G 3 = 2 

^4 “ q 3 q 2 

(2. 10c) g 4 = 2 

- q 3 q 4 q 2 

A A A 


q l q 2 q 3 _ 


q 2 -q, q 4 


Define 


enough 


( 2 . 6 ) 


(2.7) 


( 2 . 8 ) 


(2.9) 


(2. 10b) 


(2. lOd) 
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( 2 . 11 ) 


h, = G.v 

-j J-r, I 


and 


"q.i - [-1 1 -2 I -3 1 ^ 4 ] 


then using (2.12), (2.7) can be written as 


A(q)v r ^ = A(q)v r l - H q jSq 


Finally from (2.1), (2.2), (2.3) and (2.13) we obtain 


*1 = -b, i * H e, 1 S = * H e, 1-i - A( 3>* r , i * H q, 1 S 3 

Since v A(q) and v are error-free, it is clear that 
d , 1 — r , 1 


2b. i * A( S>2 r , 1 


therefore (2. 14) can be written as 


*1 = H e, i 5 £ + H q, i 5 2 + El 


where 


n* = H 
-1 


.u, 

e, i— i 


Note that (2. 16) can be written as 



T T T 

The propagation of the vector [Sq | Se ] (where T denotes the 
time can be expressed by the linear equation [8] 


Sq ' 





Se 

= 

F 


Se 

<5p _ 




_ S E . 


( 2 . 12 ) 


(2. 13) 


(2. 14) 


(2. 15) 


(2. 16) 


(2. 17) 


(2. 18) 


transpose) in 


(2. 19) 
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where 5p contains additional states necessary to express (2. 19) as a linear 
equation driven by a white noise vector w. For compatibility with (2.19), 

(2.18) is extended to include Sp as follows 

+ n 1 (2.20) 

The set (2.19) and (2.20) can be written as 

x = F x + w (2. 21a) 

Xi = H i^ + E* (2.21b) 

where H =[H |H 1 0 ] . The latter equations can be used in an EKF to compute 

x^, the estimate of x at time t^. 

T T T T A 

Let X =[q |e |p ] then according to the EKF algorithm, X^(-), the a-priori 

estimate of X at time t^ is used to calculate which is needed to obtain the 

a-posteriori estimate x^(+). The latter is then used to update the entire state 

estimate as follows 

-i (+) = *i (_) + *i ( + ) (2.22) 

Using (2.21a), X^( + ) is propagated in time to become X^ + j(-), the a-priori 
entire state estimate at time t^ + ^. The dynamics matrix for the propagation of 
X^( + ) is denoted by A (see (2.23a). The covariance which is needed for computing 
the Kalman gain necessary for evaluating x^(+), is computed according to the 
ordinary Kalman filter algorithm. To sum it up, the full EKF algorithm is as 
follows 

Between measurements 

Solve from t. to t. „ 

i l + l 

X = A[X(t), t]X (2.23a) 

P = A [X( t ) , t ] P + P A T [X(t),t] + Q(t) (2.23b) 


y. = [H . H , 
-l q, i e, i 1 


0 ] 


5q 

Se 

Sp 
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with the initial conditions X =X, ( + ), P=P. (+). The solutions at t, , are 

o i o i i + 1 

denoted by *i+l ( - } anci + respectively. Q(t) Is the spectral density matrix 
[9] of w. 

Across measurements 


W* ! - K iH.iy 1+ i 


(2. 23c) 


(2.23d) 


*i+l (+) -W’ + w +) 


W*> - [■ - W->[' - ^i+i B i+i] T * 


(2. 23e) 


(2. 23f ) 


where R is the covariance of n* .. 
i+1 -l+l 

Compensation 

In computing (2.23a) and (2.23b) we need to use the gyro output vector w 
which contains errors. Those errors are estimated as a part of p. Before each 
time (2.23a,b) are used, the errors have to be appropriately compensated using 
their estimate. Similarly, v^ which is used in (2.1) to obtain y^, contains 
erros which constitute e. Before each time y is computed (for use in (2.23d)), 
the errors in v^ m . have to be compensated using their estimate. 


III. THE ROLE OF QUATERNION NORMALIZATION 

The state measurement update given in (2.23e) can be written in an explicit 
form as follows: 


(3. 1) 


q( + ) 


q(-) 


5q( + ) 

e( + ) 

= 

e(-) 

+ 

5e( + ) 

p( + ) 


p(-) 


Sp( + ) 

- 

i+1 

- 

i + 1 

L 


i+1 
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Unless convergence has been attained, q^ + ^(+) is n °t necessarily normal even if 
q (-) is. We know, however, that the quaternion which the algorithm is trying 
to estimate is necessarily normal. We can then enforce normalization on q.^( + ) 
with the hope that the enforcement of this quality of the correct quaternion 
will direct the estimated quaternion in the right track and will enhance its 
convergence. Indeed, it was found in the past [7,8] that normalization is 
helpful. In particular, it was found that when the attitude varies very slowly 
between measurements, normalization, although not necessary, resulted in a 
faster convergence; however, when the attitude changed rapidly between 
measurements either filter tuning or normalization were necessary to avoid 
divergence. The use of normalization is superior to tuning because, first, 
tuning involves a tedious trial and error process, second, tuning is not a 
robust solution, and third, with quaternion normalization the final attitude 
estimate is closer to the correct quaternion. 

Four normalization schemes are discussed next. 


IV. BRUTE FORCE (BF) NORMALIZATION 

The BF normalization is performed as follows [7], After X^ + ^(+) has been 
computed in (2.23e) the quaternion part of the state (see (3.1)] is normalized 
as follows 

q^jOO - q ltl <*V II q 1 + l (,) 11 ,41) 


and then, the normal quaternion, q* + ^(+), i- s usec * to re-form X^ + j(+) as follows 



q* ( + ) 
e ( + ) 

p( + ) 


(4.2) 
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This straightforward mode of normalization constitutes an outside interference 
in the EKF algorithm which has to be accounted for in order to avoid filter 
divergence. It was shown [7] that the normalization operation of (4.1) is 
tantamount to the following computation of ^* + j( + ) 


■ [aui*-’ * 5 3i*i (+l ] - («-s) 

Therefore, while the EKF algorithm presented In Section II assumes that the 
a-priori quaternion estimate is updated according to (2.23e) as follows 


2 i+ l (+) = ^i + l (_) + 5 2l+i (+) (4.4) 

in reality, due to the normalization, it is updated according to (4.3). The 

difference is then in the term '^ i+1 (~^ i+1 5 g i+1 ( + ) • Because of this residual 

term, the full reset implied by (4.4) does not hold anymore. Therefore, 

following the logic of the EKF algorithm, the residual term, 

A X 

"^i + l + + 1 ^ has to ke propagated In time. It was shown [7] that this 

mode of normalization does not affect the covariance computation of the EKF; 
therefore, only the state computation has to be modified. In view of the 

normalization operation of (4.1), the following changes have to be made in the 
EKF algorithm presented in Section II. Between measurements, In addition to the 
computation of and P^ + j(“), compute also as follows. Solve from 

t<i to + ^ the differential equation 

<5q - F [X(t), t]5q (4.5) 

where F q is the 1,1 submatrix of F, with the Initial condition 

= 3* (“)qJ(~05q. ( + ) and denote the solution at t, , by <5q. , Then form 

— o il— I i+l—i+1 

Xj +1 (-) = [«§i +1 (-).0 T ,0 T ] (4.6) 

and change (2.23d) to read 
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( 4 . 7 ) 


The BF normalization algorithm has all the expected advantages mentioned in 
Section III; however, it is not elegant in the sense that the normalization 
constitutes an outside interference in the EKF algorithm which has to be 
compensated. This compensation adds a certain complication to the algorithm 
presented in Section II. Therefore we propose the following QPM normalization 
scheme. 


V. QUATERNION PSEUDO -MEASUREMENT (QPM) NORMALIZATION 

According to this algorithm the updated quaternion § i+1 ( + ) is used to form 
a pseudo-measurement as follows 


y n , i + i 


= 5 i+ i (+) /ll 3 i+ i (+) II 


( 5 . 1 ) 


It is then assumed that the quaternion is measured by an imaginary device, say 

"quaternion-meter", and the output of this device is y n> 1+1 Plus a small white 

measurement error. Following this rationale a measurement update is performed 

which is based on the quaternion measurement. To accomplish that we realize from 

(3 1) that y . , is related to the state vector as follows 
-n, x+1 


n , i+1 ^n, i+l-i+1 + n n, i+1 


( 5 . 2 ) 


where 



r T, T. T. 

[q le |p J i+1 


( 5 . 3 ) 


1 

0 

0 

0 



0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 




and R ... the covariance of n^ is the diagonal matrix 
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, 2 2 2 2 , 
= diag[r , r , r , r ] 


(5.5) 


R n, i+1 

and where r is a small number. The a-priori state estimate for this measurement 

update is, of course, X i+ j(+). Note that the output of this update is the full 

state vector and not just the estimate of, x, the difference between and 

its estimate X^ + j(+). This pseudo-measurement update is performed after the 

computation in (2.23f) has been carried out. The pseudo-measurement update 

algorithm is as follows. K is computed according to (2.23c) where H, and 

n, l+i i+1 

P i + l and R i+1 are replaced by i+J , P i+J (+) and R r i+J respectively. The 

state update is then re-computed as follows 


=W*> = W + > 


K n, i + l ^-n, i+l 


- H n,i+A+l ( + )] 


and p i+1 ( + ) is re-computed according to (2.23f) where K 1+1> H i+1 and P (-) are 
replaced by ^ + and P^ + j(+) respectively. The new estimate and its 
covariance are then propagated in time as before. 

The QPM normalization performs quite well and achieves the expected 
benefits of quaternion normalization provided r is well tuned. If this is not 
the case, the attitude estimate may reach a wrong value, and if the attitude 
changes between vector measurements, it may even diverge. The reason for this is 
described next. 

For the normalization to be effective one is tempted to choose a small r in 
which case the filter practically replaces the stored quaternion estimate by the 
normalized quaternion. However, the small "measurement noise", r, reduces the 
variance of the quaternion estimation error considerably. Therefore, the filter 
assigns a very high credibility to the normalized quaternion estimate even 
though it is not yet the correct quaternion. Consequently, the filter does not 
allow new vector measurements to alter the quaternion estimate and the latter is 
stuck on a wrong value. If the quaternion changes now due to attitude change 
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then the quaternion estimate diverges. In order to avoid this phenomenon one has 
to tune the value of r which constitutes an additional design burden. Therefore 
although, unlike the BF normalization, the QPM normalization blends naturally 
into the EKF algorithm, the required tuning constitutes a considerable 
disadvantage. To alleviate this problem we proposed the following magnitude 
pseudo-measurement normalization scheme. 


VI. MAGNITUDE PSEUDO-MEASUREMENT (MPM) NORMALIZATION 

Unlike the previous scheme, where we assumed that we "measured the 

normalized quaternion, here we assume that we "measure" the square of the 
quaternion Euclidean norm whose magnitude is assumed to be 1. This imaginary 
"norm meter" yields the reading z where 


z n, i+1 * 1 * V n, i + 1 


( 6 . 1 ) 


and where v is assumed to be a white measurement noise whose variance is r. 

n, i + 1 

Note that the "measured" quantity is a non-linear function of the quaternion 
components; therefore, we compute the effective measurement, y n> i+1 » as 


y n, i+1 = z n, i+1 ' ( 4M 1,1*1 * * ^*>3.1+1 * ,(+) 4.i+l] 


Using (6.1) and (2.4), (6.2) can be written as 

y n, i+1 1 ” I ( q j.i + l 6q j,i+l) + V n, i+1 

j=l 

Neglecting products of 5q j n+ ^, (6.3) can be written as 

y n, 1 + 1 “ 1 ' [ I q j. 1.1 ' 2 I q j. 1+I aq j. l»l] * V n. i+1 


j = l 


j = l 


( 6 . 2 ) 


(6.3) 


(6.4) 


and since 
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1 


(6.5) 


£ q j, i+1 
J=1 


then (6.4) can be written as 


n, i + 1 


[ 2q l, i+1 2q 2, i+1 2q 3, i+1 2q 4, i+l] 


Sq 

5q. 


1, i + 1 


+ v 


n, i+1 


( 6 . 6 ) 


2, i+1 


*>3, i + 1 
*>4.1+1 


Since q^ j-1,2,3,4 is unknown, we follow the common practice of replacing 

the quaternion components by their estimate, thus 


n, i + 1 


= 2 qM i,i+i 


2q( + ) 


2, i+1 


2q( + ) 3, i+1 


2q(+) 4, l+l] 


*> 1 , 1+1 
*> 2 , 1+1 
aq 3, i + 1 


+ V 


n, i + 1 


<5q 


4, i+1 


(6.7) 

The latter is the measurement equation which is used to perform a magnitude 
pseudo-measurement normalization update. The sequence of operations is similar 
to that performed when the QPM normalization update is carried out (see the 
preceding section). The only difference is that now 


y n,i+l 


-[«<*> 


2 

j. i+1 


J=1 


(6.8) 


H= [25W 1 >m | 2 q '+>2.1+l| 2 q(+) 3.1+l 


2q(+) 4,i+l] 


(6.9) 


and 


R . = r 

n, i + l 


( 6 . 10 ) 


We realize that the fact that r is very small does not imply that the 
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measurement of q is precise. It only implies that the measurement of II q II is 
precise. Therefore, now the variances of the quaternion states do not increase 
to a value close to r and thus the estimates of the quaternion components do not 
cling to wrong values and stay there like they do when the preceding QPM 
normalization is applied with a small r. 

VII. TEST RESULTS 

The algorithms presented in this paper were and still are being tested now. 
In these tests the EKF is applied to simulated as well as to real Earth 
Radiation Budget Satellite (ERBS) data. Partial results are presented as 

follows. 

Quaternion normalization speeds up the convergence of the additive EKF when 
used to estimate spacecraft attitude from vector measurements. Moreover, if the 
attitude changes considerably between vector measurements, quaternion 

normalization replaces filter tuning which is necessary to avoid divergence. In 
the latter case, quaternion normalization also reduces the final attitude 
estimation error. 

In Table 7. 1 we see the final attitude estimation error when the EKF is 
applied to simulated ERBS data. The initial attitude error is 30° and the value 

-5 

of r used in the QPM and MPM algorithms is 10 . 


Table 7.1: Final Attitude Error in Degrees at 100 sec, r-10 



Normalization Algorithm 


Without 

Normalization 

BP 

QPM 

MPM 

Yaw 

.0048 

.0074 

.0057 

. 0069 

Roll 

.0022 

-. 0002 

.0019 

. 0039 

Pitch 

.0170 

.0060 

-. 0009 

-. 0033 

RMS 

0. 0178 

0.0095 

0.0061 

0. 0086 
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Note that the BF algorithm implies no measurement update therefore no r is used 
in this run. We turn to Table 7.2 to see the advantage of the MPM over the QPM 
algorithm. We realize that while for r=10~ 5 both algorithms exhibit identical 
accuracy, the QPM algorithm fails when r=10 _11 . The reason for this difference 
was mentioned at the end of Section V. 

Table 7.2: Final Attitude Error in Degrees at 100 sec, r=10 



Normalization Algorithm 


QPM 

MPM 

Yaw 

3. 2387 

.0045 

Roll 

10. 3660 

.0083 

Pitch 

-0. 7451 

.0127 

-..I. 


VIII. SUGGESTED FUTURE RESEARCH 

Although the MPM normalization performed satisfactorily we suggest to 
investigate an algorithm of implied normalization which does not really use 
normalization. This algorithm is presented next. 

In Section II we presented the development of the additive EKF for 
quaternion estimation. In that development we derived the linearized 
relationship between the vector measurement error and the quaternion estimation 
error which are summarized in (2.16). To meet this end we differentiated the 
matrix A(q) given in (2.9). The differentials were partial differentials with 
respect to the elements of q. As a result of the differentiations we obtained 
the matrices G^, j=l,2,3,4 which are listed in (2.10). 

When q is indeed of unit length, A(q) is an orthonormal matrix; that is, 
its columns (rows) are orthogonal to one another and are of unit length. If, 
however, q is not of unit length, then the columns (rows) of A(q) are still 
orthogonal to one another, but their length is not a unit anymore. It was proven 
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in [10] that the matrix A*(q) computed as follows 


A*(q) = 


II" q 


-2 A(q) 


( 8 . 1 ) 


is not only orthonormal, but it is also the "closest' 1 orthonormal matrix to 
A(q); that is, of all possible orthonormal matrices, the distance between A*(q) 
and A(q) is the smallest where by distance we mean the Euclidean norm of the 
difference matrix A*(q)-A(q). It can be argued that if we use A*(q) rather than 
A(q) , we practically enforce normalization. This is so because normalizing q 
first and then using the normalized quaternion to compute A(q) is identical to 
the computation of A*(q) as given in (8.1). The partial differentiation of (8.1) 
with respect to the quaternion components yields 


5A*(q) 


r j - r j ( 3> * aq 


J 


2q. 


II q II 


A(q) + 


4 


II q II 


2 G J ( 2> 


( 8 . 2 ) 


where Gj(q) is given in (2.10). The final algorithm is as given in Section II 

with r. replacing G. in (2.11). We call this normalization scheme the linearized 
J J 

orthogonal i zed matrix (LOM) algorithm. 

Finally, in the future we intend to apply all the normalization schemes 
discussed here to real ERBS data. 


IX. CONCLUSIONS 

It was found again that quaternion normalization in the additive EKF for 
attitude determination from vector measurement has the following advantages. If 
the attitude changes slowly, normalization speeds up estimation convergence. If 
attitude changes rapidly between measurements and no normalization is applied 
then filter tuning has to be used in order to avoid divergence. However, if 
normalization is applied, convergence is achieved without filter tuning. 
Moreover, the final attitude estimation error is smaller. There is then a clear 
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advantage to quaternion normalization. Three quaternion normalization algorithms 
were tested. The conclusions with regard to the use of each one of them is 
listed next. 

■ The brute force (BF) normalization algorithm works well and exhibits the 
normalization benefits described before. 

■ The quaternion pseudo-measurement (QPM) algorithm performs well only after 
tuning. 

■ The magnitude pseudo-measurement (MPM) algorithm performs well and needs no 
tuning. 

Finally, we suggest the investigation of the linearized orthogonal matrix 
(LOM) normalization whose development was presented in Section VII. All the 
normalization schemes will be tested on real ERBS data. 
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